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Objective 
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This report describes work performed under co-operative agreement NCC2-767 with 
NASA Ames Research Center, during the period from June 1993 to January 1995. The 
objective of this research was two-fold. The first objective was to complete the three- 
dimensional unsteady calculations of the flow through a new transonic turbine and study 
the effects of secondary flows due to the hub and casing, tip clearance vortices and the 
inherent three-dimensional mixing of the flow. It should be noted that this turbine was 
and is still in the design phase and the results of the calculations have formed an inte- 
gral part of the design process. The second objective of this proposal was to evaluate the 
capability of rotor-stator interaction codes to calculate tonal acoustics. 


Motivation 

The ultimate motivation behind this proposed research is to be able to simulate a com- 
plete propulsion system. In order to do this, certain key CFD methodologies such as 
turbomachinery flow solvers have to evolve to a level of maturity so as to be used with 
confidence. An important criterion for a numerical code to be a design tool for turbo- 
machinery applications, is its ability to predict accurately, the secondary flow features, 
losses and other more sensitive flow quantities such as tonal acoustics. 


Background 

Computational fluid dynamics (CFD) is playing an increasingly important role in the de- 
sign of various propulsion components. Considerable progress has already been made in 
using CFD in the design of turbopumps and impellers. Turbomachinery (rotor-stator) 
flow solvers have been developed at NASA/ARC which incorporate some of the most 
modern, high-order upwind-biased schemes for the solution of the thin-layer Navier- 
Stokes equations (Ref. 1-2). The family of rotor-stator interaction codes (which includes 
ROTOR- 1, ROTOR-2-4 and STAGE-2 codes) are currently in use in several industrial 
and government organizations. So far, the unsteady rotor codes have proven quite capa- 
ble of predicting pressure variations on the surface of the airfoils. They have also demon- 
strated their capability of predicting more “sensitive flow variables such as the total 
pressure losses in the flow field. In addition, two-dimensional versions have been used in 
the actual design of turbomachinery (Ref. 3). However, the success of CFD in the de- 
sign of turbomachinery has largely been in predicting the two-dimensional time-averaged 
and unsteady pressure loads on airfoil surfaces. However, there are important three- 
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dimensional effects in the flow associated with turbomachines such as secondary flows 
due to the hub and casing, tip clearance vortices and the inherent three-dimensional 
mixing of the flow which require additional detailed analysis. Additionally, the flow in 
the tip clearance region is not well understood. There can be considerable flow turning 
and temperature variation in this region that can affect the overall performance of the 
turbomachine. An important aspect of the flow in a turbomachine is the impact on the 
overall losses due to secondary flow features. Hence the timely completion of the three- 
dimensional unsteady flow in a new transonic turbine in the design stage would make 
rotor-stator interaction codes developed at NASA/ARC a viable tool for turbomachinery 
design. 

The second aspect of this proposal was the calculation of tonal acoustics in turboma- 
chines. In Ref. 4-5, two-dimensional unsteady rotor-stator interaction calculations were 
performed to study the plurality of spinning modes that are present in such an inter- 
action. The propagation of these modes in the upstream and downstream regions was 
analyzed and compared with numerical results. It was found that the numerically calcu- 
lated tonal acoustics could be affected by the type of numerical boundary conditions em- 
ployed at the inlet and exit of the computational boundaries and the grid spacing in the 
upstream and downstream regions. Results in the form of surface pressure amplitudes 
and the spectra of turbine tones and their far field behavior were presented. The mode- 
content” for different harmonics of blade-passage frequency was shown to conform with 
that predicted by a kinematical analysis. It was however assessed that a similar three- 
dimensional calculation would require a highly accurate algorithm since relying on very 
fine grids would be impractical. Also, three-dimensional non-reflective boundary condi- 
tions would have to be developed. It was with the above in mind, that the development 
of a new high-order accurate multi-zone Navier-Stokes code was initiated. This code is 
based largely on the ideas presented in Ref. 6. Figure 1 shows a comparision of the re- 
sults obtained by the new code with that obtained from ROTOR-2. These preliminary 
results look quite promising and have indicated a further study in the development of 
the new method. 


Achievements 

Flow predictions in an advanced transonic turbine was completed in a timely fash- 
ion. The task for the calculation was given at the same time as fabrication was initiated. 
The numerical results were obtained well in advance of the first experimental runs and 
have already played an itegral part in determining the placement of the probes and have 
also facilitated in the understanding of the experimental results. This exercise has made 
three-dimensional rotor-stator interaction codes a viable tool in the design process. 

A detailed study of the capability of the two-dimensional rotor-stator codes in com- 
puting tonal acoustics was completed. In anticipation of extending this capability for 
three-dimensional predictions, development of a new high-order-accurate flow solver was 
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initiated. Preliminary results appear to be promising. 
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the linear behavior <il the How in (In- fur field Two I s t ».-<! of noiircltcrlive l.onndurv ion. ill ions i,. s ied 1 1„. |j rst j s a dif 

ferenli.il oni-dimension onndar) clvriv.d f.om the far liel.l a-cmstiral behavior of fl„- (low ami similar l<> that 

developed hv llayli*. and Turk-'l t[0). |I0) and |l l|) I he second was an approxi two dimensional unsteady boundary 

condition developed m |I2|. Hits second approa.li (using different nonreflecl ive boundary condii.ons) will he denoted the 
"short- grid’ approach l or the purpose of comparison, short grid calculations with reflective boundary conditions .ire also 
performed. 

This study focusses on an examination of the modes present in a mlor-M alor interaction for a ;i stalnr/Trolor rase 
usmg n 4 nvior-Siokt'^ solution proofin' Results m (lie form of (url'ino tom* sperJra mol axi:il variaiinn of amplitude* of 
sonii» inodes ;iit pr»-<«'ni t c! 

KINEMATH'AI. ANALYSIS 01 IN I'! ill AT Tit IN TON K NOISE 

Reflective ho.indary conditions, such as a lixed exit static pressure condition, can lie made to behave essentially as 
noiireflcrtive boundary conditions will) the use of appropriate grids in (lie far field region When the grid cell sizes in the 
far-lieUI region arc of the order of the wavelength of the mode to be attenuated, the energy associated wit), that mode de- 
cays rapidly because of numerical dissipation. Hence, in order lo atlcnuale relleclions at. the computational boundaries 

when leflective boundary conditions are used, an estimate maxim wavelength in the pressure field is required. 

here fore a kinematic,-, I analysis of the I'm. tier modes was earned mil This analysis, when coupled will, a linearized anal- 
ysis of (he flow gives important information about the relative magi, it tides (of the diflWcnl Fourier modes) as a function of 
distance. I he method was first used in (fij for the analysis of a single stage and can he generalized for multiple stages The 
assumption made here is that the tone generating mechanisms occur at multiples of Made passing frequency It should he 

noted that contributions from wake shedding and other ae.odynamic noise sources are ignored The results seen. 

to indicate that these secondary conlrihi, lions are small for the geometry I (low conditions chosen. 

I he billowing analysis is limited to two-dimensional How hut can he extended lo three dimensions in a straightforward 
fashion Consider a single-stage turbine, where S is the number of stator airfoils and It (lie number of rotor airfoils. The 
composite pitch is the transverse distance over which the Mow is periodic and is denoted by I The distance between the 
stator airfoils is I/S and the distance between rotor airfoils is I/ll. The velocity of the rotor is denoted by V„ The pres- 
sure at any axial plane in the flow field, is assumed lo he periodic in time will, period equal lo the minimum i imr required 
for the rotor-stator geometry to repeat In the case of a single-stage configuration with identical, equispaced rotor airloils 
Hie period in a stationary frame of reference is given by l/\V„\l{. This is also know,, as the blade passing lime The pres-’ 
sure variation in any particular axial plane is given by 

= E E " «»(y(»»y-»i«|FR|0 + * ) (1) 

where u is the harmonic of blade passing frequency and m corresponds to the spatial I, anno, lie in „ Assuming that ev- 
ery stator airloit the stator row js equally spaced, the pressure variation ,s assumed lo satisfy the shift conditio,,, which 
slatis t i.it Pi!/, t) — p(y - s , l — which is valid only when frequencies that ate non-coiimieusurafc with the blade 

passing frequencies are absent. This then yields 


m = n R(sgn{ V n )) ± kS (2) 

whore b - .... -1,0, 1, ... is the spatial harmonic of the disturbance produced Ivy the stators 
Propagating Modes 

Analytical solutions representing the unsteady flow in the far field can also he derived, as given in (8) It. is assumed 
that in the far upstream and downstream regions, the unsteady flow is a linear perturbation or a steady uniform How. (The 
underlying steady uniform How in the upstream region is different from that in the downstream region) A linearized solu- 
tion of the pressure field indicates that a Fourier mode propagates, if 


’>K Mn . M v 

7 r ~ ~ -= -f — . ~ > 1 

'» v/1 - M] v/T 


( 3 ) 


where M/i is 1^1, Af, is the axial Mach number of the underlying How and Ah, is the transverse Marl, number of the un- 
derlying flow. I lie axial wavelength of the propagating mode is denoted by A, and is given by 


/(I ~ A I ; ) 

A/,(>i ft Mr, + A/„ in) ± ^'HVn)^/(nliM n + M tJ m ) 2 - (1 - A IJjnd 


('!) 


where m is given by Eq. (2). For very low values of axial and transverse Mach numbers. Af r . and AT respectively the 
maximum wavelength A r cat, he approximated by substituting 0 for A/, and A/„ in Eq. (•!). For low values of it cnil 
1)0 seen Mint Tor the lovver liar monies. 


IM. 


i 

>i(l. Mu 
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For a 3-sta lor/ Trot or case that is presented hero, (A,] IBor = I/IVIM,,) since the fundamental and the first harmonic do 
not propagate under the present assumptions of AT, Af„ « I When a long grid is used lo dissipate the propagating modes 
the far held region, the grnl sparing' near the exit and inlet boundaries are chosen to he about [A,) r /2 it should be 
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noted that there is a possibility of large axial wave lengths when u II M fi a; \J ( I — A/ 2 ),n, but (bis would occur only for 
higher harmonics under I lie present assumptions of low Mach numbers. 

Decaying Modes 

A fourier mode will decay with increasing axial distance from I In- rotor stator pair if the inequality sign in Eq, (3) is 
reversed. The amplitude of the mode would vary as 


= (y jtittH * 1 - - \" RMn h ,,,A M 


A decaying mode decays exponentially. It is understood that for z < 0 the positive sign in Eq, (0) is used and for z > I) 
the negative sign is used. The variation of the natural logarithm of the amplitude (log(«„,„)) is linear with slope 

GEOMETRY AND GRID SYS TEM 

The airfoil geometry used is that given in |2). The two-dimensional compulations ol tins study were performed using 
the experimental airfoil cross-sections at mtdspan. These cross-sections, along with a schematic of tin' computational grid 
arc shown in Figure !. A system of patched and overlaid grids is used lo discretize the How region of interest. The inner 
grids are O-grids that were generated using an elliptic grid generator. The outer li-grids were generated algebraically 1 lie 
experiment consisted of 22 stator airfoils and 28 rotor aitfoils. To model the experimental setup the How over at least 25 
airfoils (1 1 stator airfoils and M rotor airfoils) would have lo he calculated. This would require exresslvc computational 
resources It was therefore decided to solve a smaller problem by using rescaling strategies as shown in (I3| in order to re- 
duce the airfoil count. The number of stator airfoils was changed from 22 lo 21 and the size of each individual stator airfoil 
was enlarged by a factor of 22/21 This rescaling allows a 3-slalnr/ l-rotor computation wherein periodicity is imposed on 
the flow over 3 stator airfoils and 4 rotor airfoils. Clumping the airfoil count does change (lie nature of the tonal acoustics 
in the flow Held because the mode content of the propagating modes depends upon the airfoil count as indicated |,y Eqs. 

(2 7). However, since the objective of this preliminary investigation is lo evaluate the capability of rotor-stator interac- 
tion codes ; u calculating tonal acoustics, the rescaled rotor slatm geometry was used. It should he mentioned that figure 
I shows only a schematic of the grid The ai tinl number of grid points is lunch latger and I In: spacing between the grid 
points is mmli smaller. More details about I lie grid system can hi- found in [l.lj 

NFMI'lt If ’A I METHOD 

The unsteady, thin layer, Navier Stokes equations are solved using an upwind- hia-ad Imilc-dill.-rciicc algorithm. The 
kiln-main viscosity was calculated using Sol herfands law and tin luilniteiil »ddv vi-ao-uly was calculated using the 
Hahlwin l.oinax model. I he method is Ihiid-onb-i accurate in spa, e and second older-accurate in time. At each lime step, 
several Nr uoai delations aie performed, s,. that tin- hilly implicit finite difference equations are solved. Additional details 
regarding tin seln-me c an be found in |I3], 

UODNDAHY CONDITIONS 

Tin I •oinnlai y conditions required win'll using multiple zones .an In- broadly classified into Iwo types. I lie first type 
Consists of tin- zonal conditions which arc implemented at the inlcifa.es of I in- computational grids and I lie second type 
consists of tin- natural boundary conditions imposed on the surface- and tin- ontei boundaries of tin* computational grid. 

The treatment of the zonal boundaries ran lut found in [I l|. The natural boundary conditions used in this study arc dis- 
cussed below. 

Airfo il Snrfa cc Don ndary 

1 lie- I iciii it i l;i i v enuditions on the airl'oii surfaces are the "no slip" rendition and adiabatic wall conditions. It should be 
noted that in the case of the rotor airfoil, ‘no-slip" implies zero relative velocity at the surface of I lie airfoil. In addition lo 
the ‘no slip" condition, the. derivative of pressure in the direction normal -to the wall surface is set In zero. 

Exit Houndary 

One reflective and two radiative boundary conditions were studied, for the reflective boundary condition (see e.g. 

[13]) the exit pressure was specified and three quantities are extrapolated from the interior. The three quantities are the 
Uiemami variable R\ = u + 2 c/(y - 1), the entropy .S' = | >//>' and ti the transverse velocity. This type of boundary condi- 
tion reflects the pressure waves that reach t lie boundary back into the system. 'Two types of radiating boundary conditions 
were also implemented. The first was a one-dimensional boundary condition formulated by llayliss and lurkel ([If] and 
also [18)) ft is assumed that at the downstream boundary, the flow is linear. Two-dimensional boundary conditions as pre- 
sented in (12) were also implemented. As in the previous boundary condition, liic flow at the exit is assumed to have small 
perturbations and lienee linearizahle about an underlying mean flow. Implementation of tliis boundary condition requires 
a knowledge of the underlying exit flow variables, p^., u, 0 . i, N , and Pco ■ The first three quantities are time lagged whereas 
the exit pressure. p M , is kept constant. 

Inlet llcnudary 

One reflective and two radiative boundary conditions were also studied for tin' inlet. The first was lire reflective 
boundary condition procedure wherein three quantities have to In' specified. 'The three chosen are the Ricmann invariant 
/?!=« + 2r/(q - I) = ii(-cc) + 2c ( _ c<: ,/(-> - 1) the total pressure ;i co <a/ = P(_v,)(l + (7 - l/2)A/ ( J _ rc) ) ^ and the inlet 
flow angle, which in this case is equivalent to v,„i rt - 0 The fourth quantity needed to update the points on this boundary 
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ABSTRACT 


Turbochargers, with maximum running speed in excess of 200,000 rpm, 
are perhaps the fastest running turbomachinery in today's Automotive 
industry. Turbocharger noise problem came into prominence with the 
application of turbochargers on passenger cars in early 1970. During 
the early period of its application, a "pure tone" noise synchronous 
with the turbocharger speed, was the dominant objectionable noise and it 
was controlled by reducing the imbalance level. Since then, additional 
types of turbocharger noise, excited by mechanical vibration and ait 
flow have been encountered . The trend towards ) ightet and more compact 


vehicles with reduced engine space has increased vehicle's response to 
turbocha t ger generated vibration and pressure pulsations that result in 
vehicle noise. This presents an ever-increasing challenge to the 
turbocharger manufacturers tor producing quiet turbochargers. 'I his 
paper presents a brief description of the different types of 
turbocharger noise, their generating mechanisms and control. In most 
cases, the noise reduction is achieved by design improvement and better 


control of manufacturing process. This consists of changing the rotor- 
bearing system design, improving the balance, process to correct 
extremely small amount of imbalance and tightening the casting and 
machining tolerances for more symmetrical compressor and turbine wheels. 


INTRODUCTION 


Turbochargers are very high speed turbomachinery, used primarily on 
automotive vehicles to increase the power output of a given size, of 
engine. Till the end of last decade, turbocharger noise was addros ,u 1 
on a "fire fighting" basis. A prior assessment of potential noise 
problem was not done during the design and development phase which 
resulted in high cost for noise control . This has now started to change 
and potential noise problems are being addressed from the start of 
design. This paper will describe the various noise generating 
mechanisms and the steps taken during the development of a new 
turbocharger to control the noise generation. The development of this 
turbocharger is nearing completion and the information presented is 
based on development testing. The next step in the cost effective noise 
control is for the turbocharger manufacturer to provide turbocharger 
related input to the vehicle manufacturer at the start of a new 
engine/vehicle program and stay involved till the design is complete. 
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In this study, the unsteady, thindayer Navter-Stokes equation^ solved ‘ ^g^system 

SflMilllll 


1. INTRODUCTION 

Tm FLUID flow within a turbomachine is inherently unsteady. There are several 
^chankms that cause the unsteadiness, and some of these are, (a) the relative motion 
between the rotors and stators (which is also called the invisctd or potential effect), 
(b) the interaction of the downstream airfoils with the wakes generated by th 
upstream airfoils, and (c) the shedding of vortices at blunt trailing edges^In general as 
the axial gap between the stator and rotor airfoils decreases, the magnitude ot the 
unsteady interactions increases. These interactions can even become strongly coupled. 
Hence, to study unsteady processes involved within a turbomachme, it may be essential 

tn t-rpflt the rotor and stator airfoils as a single system. , 

Pioneering work in predicting inviscid rotor-stator interaction was conducted by 
Erdos et al .\ s far back as 1977. However, the subject has only very recently besom 
the focus of increasing attention due to the considerable increase in computationa 
resources Lewis et al (1987) solved the quasi three-dimensional mviscid equations, 
Jorgenson & Chima (1988, 1989) used the explicit Runge-Kutta metho ^ “ S °^® ® 

quasi-three-dimensional Euler and thin-layer Navier-Stokes e q ua b ons k ™ re 
dimensional Deriodic calculations have also been presented by Koya & Kotake (19® )• 
Gibeling et al. (1986) presents results for the flow in a compressor stage obtained using 
^shearing grid technique. Here, a single grid is used to discretize the flow domain and 
t allowed to shear in order to allow relative motion between the rotor and stator 
airfoils. The data from the sheared grid are interpolated onto an undistorted initial gn 
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at regular time intervals so as to limit the cumulative distortion due to grid shear 
The development of general zonal techniques and robust, accurate algorithms for 
numerical solution of Euler and Navier-Stokes equations has contributed further to 
development of rotor-stator interaction codes. Giles (1988a) has calculated t 
dimensional rotor-stator interactions using the Euler equations. In this work, a nc 
concept of a “time-inclined” computational plane is used in order to s’urmo 
difficulties encountered when the stator-rotor pitch ratios is not a ratio of two sr 
integers. Rai (1987) presented a two-dimensional calculation of rotor-stator interact 
for an axial turbine. The airfoil geometry and flow conditions were those given in Dr 
et al. (1982). More recently, Rai (1989) and Madavan et al. (1989) computed the fi 
three-dimensional flow fields for the same case. The hub, outer casing and the rotor 
were a11 deluded in the calculation. Rai (1987, 1989) and Madavan er 
( 989) solved the thin-layer Navier-Stokes equations in a time-accurate manner us 
an implicit, upwind-biased, third-order accurate method to compute the flow field. 1 
ability of their codes to predict near field flow quantities, such as the time-averat 
pressure distributions on airfoil surfaces and the pressure amplitudes and phase on "i 
surface of the airfoils, was demonstrated. In addition, the two-dimensional codes w< 
used to predict accurately the total pressure defects in wakes. These compu 
programs have also been recently used in the design process of turbomachir 
(Kangwalla et al. 1992). Rotor-stator interaction codes have also been used to calcuk 
sensitive flow quantities such as heat transfer (Rao et al. 1992a, b; Griffin 

^ C fn°r!^ aU8hey 1989 ^ and three * dim ensional unsteady hot streak migration (Dorney 
al. 1990). Recently, the codes have also been extended for multiple stage calculatio 
by Gundy-Euriet et al. (1991). More recently, Dorney (1992) performed a rigoro 
validation of a modified rotor-stator algorithm through comparisons with analytical a: 
mearized unsteady aerodynamic solutions. The present work is focussed on investig; 
mg the ability of rotor-stator codes to predict “correctly” the tonal acoustics in t 
now field due to rotor-stator interactions. 

An axial flow turbine (or compressor) produces rotating pressure patterns call, 
spinning modes that may propagate in a spiral path. In two dimensions, these spinnii 
modes propagate at a non-zero angle to the axial direction. For any particul 
harmonic of blade passing frequency, the interaction field can produce an infini 
number of spinning modes. Each of these modes rotates at a different speed. Some . 
these modes propagate, whereas others decay. The modes that are possible for son 
multiple of blade passing frequency depends upon the number of stator and rote 
airfoils. This study focusses on an examination of the modes present in a rotor-state 
interaction for both a single-stator/singie-rotor and a 3-stator/4-rotor case using 
Navier-Stokes solution procedure. 

The numerical study of tonal acoustics involves the study of the flow field in tb 
upstream and downstream regions of the interacting rotor-stator airfoils. This raise 
t e issue of the numerical boundary conditions employed at the computational ink 
and exit boundaries and the grid spacing used in the upstream and downstream region: 
raditionally, the boundary conditions employed have been reflective. Reflectiv 
boundary conditions have been used because they provide greater control on turbin 
operating conditions, such as mass flow and pressure ratio. However, reflectiv 
boundary conditions are generally inadequate for the study of tonal acoustics since the 
reflect the propagating modes back into the flow field. Rai & Madavan (1990 
performed a 1-stator/l-rotor calculation with non-reflective boundary conditions basei 
on Riemann invariants. These non-reflective boundary conditions were the same a 
those developed by Erdos et al. (1977). However it was found that the use of sucl 



rotor-stator interaction acoustics 


613 


• u. j •_ a , n ^ 0 f control over the mass flow rate through the 

boundary conditions resulted in , tried. The first approach was 

turbine. To overcome this prob em, two app ^ computational domain 

to use reflective boundary of the 

upstream and dowistr^m of A^rotOT-s^ator^ ^ ^ determine grid spacmgs at the 

wavelengths of the propagating n or ooaeating modes are numerically 

far upstream and downstream regions, so that th ^propagat^g ^ methodology 

dissipated near the boundanes. n ^ou^ ^ and Verdon (1989) is applied to 

developed by Tyler & Sofnn (19 ), modes This approach will be referred to 

determine the wavelengths of the propaga 1 g * use a s ^ ort grid, but with 

as the “long-grid” approach. 'Hie secon *PF_^ measure of control over the 

non-reflective boundary conditions, w P conditions are based on the linear 
turbine operating conditions. These oun ary fi d boundary conditions 

&Turtel f (lm 1982a, b). The second was at dffierem 
boundary condition developed by Giles (1988W- approach. For the 

•WSS-- «- 1.*- 

configurations: a single-stator/single-rotor case short-arid approach 

cases were computed using the long-gnd approach a well l as “n PP 

and the corresponding results are “-nP^fAdd.' honall lyj result « 

SmifflrWU of the different modes are presented. 

2. K1NEMATICAL PREDICTION OF INTERACTION TONE NOISE 
Reflective boundary condittons. such as a fixed exit static P^or^iom ^ * 

analysis, when coupled with a linearized analyse of ^ ho ^ * 
information about the relative magmtu es Sofrin (1970) for the 

=s5Tr 

blade passing frequency. It should Be notea t , d results se em to indicate 

- - — 

Ch0Sen ' , • t f „ llnw . is limited to two-dimensional flow. However, extension 

,0^:1— f °s‘ straightforward. Consider a single-stage turbine as shown in 



614 


A. A. RANGWALLA AND M. M. RAI 



Direction of motion 
for rotor 



(b) 



-50 Chords 
upstream 


50 Chords 
downstream 


Figure 1. (a) Rotor-stator geometry; (b) short grid; (c) long grid. 


Figure 1(a). Let 5 be the number of sator airfoils and R the number of rotor airfoils 
The composite pitch is the transverse distance over which the flow is periodic and i 
denoted by /. Hence, the distance between the stator airfoils is L/S and the distance 
between rotor airfoils is l/R. The velocity of the rotor is denoted by V R . The pressure 
at any axial plane in the flow field, is assumed to be periodic in time, with period equa 
to the minimum time required for the rotor-stator geometry to repeat. In the case of < 
single-stage configuration with identical equispaced rotor airfoils, the period in t 
stationary frame of reference is given by U\V r \R. This is also known as the blade 
passing time. The pressure variation in the axial plane can be written as 


p(y> 0=2 Pn(y) 

// = 0 


cos 


Ixnt 


\Vr\R 

l 


+ *n{y) , 


a: 


where p n {y) is the amplitude of the nth Fourier component and 0 n (y) is the phase. It 
should be noted that p„(y) and 0„(y) are periodic functions in y with period equal to 
the composite pitch /. Using this periodicity in y and appropriate trigonometric 
identities, ultimately results in 


‘x. tc 

P(XU)=2 2 Pmn 


n— 0 m = — * 


(2) 
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where 


= a r 


cos 


~2it 

_T 


(my - nR \V R \ t ) + 0, 


(3) 


where n is the harmonic of blade passing frequency and m corresponds to the spatial 
harmonic in y. The next assumption that is made is crucial to the analysis. Considering 
that all stator airfoil in the stator row are equally spaced, the pressure variation is 
assumed to satisfy the shift condition, which states that 

p(y,t)=p(y- L s ,t-~). (4) 


It should be noted that both equations (1) and (4) are valid only when frequencies that 
are non-commensurate with the blade passing frequencies are absent. Substituting 
equation (4) in equation (2) then yields 

m - nR sgn(V/?) ± kS (3) 

where k = 1, 0, 1, ... is the spatial harmonic of the disturbance produced by the 

stators. 


2.1 Propagating Mode 


Analytical solutions representing the unsteady flow in the far held can also be derived 
(Verdon 1989). It is assumed that in the far upstream and downstream regions, the 
unsteady flow is a linear perturbation of a steady uniform flow. (The underlying steady 
uniform flow in the upstream region is 'different from that in the downstream region.) 
A linearized solution of the pressure field (Verdon 1989) indicates that a Fourier mode 


propagates, if 


nR M r ( M y 
T^vi - Mi Vi ~M l x 


( 6 ) 


where M R is \V R \lc, M x is the axial Mach number of the underlying flow and M y is the 


transverse Mach number of the underlying flow, 
propagating mode is 


/Cv = 


1 


1 — Ml 


M. 


-±sgn(V R ) 

c 


The axial wave-number of the 


— a~ 


(7) 


where 


2k 

j 8 = y (nR \V R \ + cM y m), 


2 Km r- —2 

a = Vl -Ml 


and c is the sonic velocity in the underlying steady flow. The axial wavelength of the 
propagating mode is denoted by A* and is equal to 2Klk x . Using equation (7) we get 

l(l~ ml) 

Ajt ~ M x (nRM R + M y m) ± 

where m is given by equation (5). For very low values of axial and transverse Mach 
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numbers, M x and M y , respectively, the maximum wavelength, A*, can be approximate 
by substituting 0 for M x and M y in equation (8). For low values of M R , it can be see 
that for the lower harmonics, 


[^x]max D w i (- 

nRM R 

In the case of a singie-stator/singie-rotor calculation, all harmonics would have spati; 
modes that would propagate as indicated by equation (6) and hence the maximui 
possible wavelength at the lower harmonics is 


[A,] 


max 


Mr * 


whereas, for a 3-stator/4-rotor case. 


[A,] 


max 


l 

3 RM r ’ 


(10£ 


(lOt 


since for this case the fundamental and the first harmonic do not propagate under th 
present assumptions of M x ,M R «l. When a long grid is used to dissipate th 
propagating modes in the far field region, the grid spacings near the exit and inle 
boundaries are chosen to be about [A. r ] max /2. It s hould be noted that there is 
possibility of large axial wavelengths when nRM R = V(1 — M 2 t )m, but this wouid occu 
for higher harmonics under the present assumptions of low Mach numbers. 


2.2. Decaying Mode 

A Fourier mode would decay with increasing axial distance from the rotor-stator pair 
the inequality sign in equation (6) is reversed. The amplitude of the mode would var 
as 

a mn a t ±dmnX , (11 

where 


dL, = - "?)(™) 2 - 1»RMr + (12 

A decaying mode decays exponentially. It is understood that for x < 0 the positive sigi 
in equation (11) is used and for „r>0 the negative sign is used. The variation of th< 
natural logarithm of the amplitude, log (a„„), is linear with slope ±d„„. 


2.3. Multiple Stages 

The analysis is similar in the case of multiple stages. Let the number of rotor rows be 
and number of stator rows be s. Denote the number of airfoils in the rotor rows a 
R\, Rz, • • • , R r and the number of stator airfoils in stator rows as Si, S 2 , . . . , S s . The 
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composite pitch is again denoted by l, so that the distance between the airfoils in rotor 
row i is l/Rj, and the distance between the airfoils in the stator row i is l/S f . Define 


and 


R = Highest common factor (R u R 2 , • ■ ■ , R r ) 
S = Highest common factors (Si, S 2 , , S s ). 


The stator and rotor rows can be aligned arbitrarily; however, it is assumed that all the 
rotors have the same velocity, V R . It is also assumed that in each individual row, the 
airfoils are identical and equispaced. The minimum time for the rotor-stator geometry 
to repeat is again equal to H \V R \ R. Thus, the pressure variation on any axial plane is 
given by equation (1), and the shift condition is given by equation (4). R and S can now 
be considered to be the number of rotor and stator airfoils on an “equivalent” 
single-stage configuration. The rest of the analysis is similar. 


3. GEOMETRY AND GRID SYSTEM 

The airfoil geometry used is that given by Dring et al. (1982). The two-dimensional 
computations of this study were performed using the experimental airfoil cross-sections 
at midspan. These cross-sections are shown in Figure 1(a). A system of patched and 
overlaid grids is used to discretize the flow region of interest. Figure 1(b) shows a 
typical system of grids used in this study. The inner grids are O-grids and were 
generated using an elliptic grid generator. The outer H-grids were generated 
algebraically. The experiment consisted of 22 stator airfoils and 28 rotor airfoils. To 
model the experimental set-up the flow over at least 25 airfoils (11 stator airfoils and 14 
rotor airfoils) would have to be calculated. This would require excessive computational 
resources. It was therefore decided to solve a smaller problem by using rescaling 
strategies adopted by Rai (1987) and Rai & Madavan (1990) in order to reduce the 
airfoil count. 

The tonal acoustics for two different rescalings were computed. In the first case, the 
number of rotor airfoils was changed from 28 to 22. The size of each individual rotor 
airfoil was enlarged by a factor of 28/22 such that the pitch-to-chord ratio of the rotor 
was unchanged. This rescaling results in a turbine which has equal number of stator 
and rotor airfoils, thus allowing a single-stator/single-rotor calculation. (Flow periodi- 
city is imposed over one stator airfoil and one rotor airfoil.) In the second case, the 
number of stator airfoils was changed from 22 to 21. The size of each individual stator 
airfoil was enlarged by a factor 22/21. This rescaling allows a 3-stator/4-rotor 
computation wherin periodicity is imposed on the flow over three stator airfoils and 
four rotor airfoils. Changing the airfoil count does change the nature of the tonal 
acoustics in the flow field because the mode content of the propagating modes depends 
upon the airfoil count as indicated by equations (5-12). However, the objective of this 
preliminary investigation is to evaluate the capability of rotor-stator interaction codes 
to calculate tonal acoustics. Hence a numerical solution of a rescaled rotor-stator 
geometry can be used to establish numerical boundary condition and grid require- 
ments. Figure 1(b) shows a typical grid for the single-rotor/single-stator calculation. 
However, if the grid spacing at the upstream and downstream boundaries is chosen to 
attenuate reflections as described in the previous section, the grid would have to be 
lengthened as shown in Figure 1(c). It should be noted that both the long and the short 
grids are identical in the near field region. 
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3.1. Grid Density 

It should be mentioned that Figure 1 only shows a schematic of the grid. The actual 
number of grid points is much larger and the spacing between the grid points is much 
smaller. In all the calculations, each inner O-grid had 151 points along the airfoil 
surfaces and 41 points in the wall normal direction for a total number of 6,191 grid 
points in each O-grid. Each H-grid had 71 grid points in the y-direction and 90 to 141 
grid points in the x-direction (90 points for the short grid case and 141 points for the 
long grid case). This results in each short H-grid having 6,390 grid points and each long 
H-grid having 10,011 grid points. The total number of grid points used for the 
3-stator/4-rotor short grid computation was 7 X (6,191 + 6,390) = 88,067 whereas for 
the long grid computation, the total number of grid points was 7 X (6,191 + 10,011) = 
133,414. Rai and Madavan (1990) give more details about the grid system. 


4. NUMERICAL METHOD 

The unsteady, thin-layer, Navier-Stokes equations are solved using an upwind-biased 
finite-difference algorithm. The method is third-order-accurate in space and second- 
order-accurate in time. At each time step, several Newton iterations are performed, so 
that the fully implicit finite-difference equations are solved. Additional details 
regarding the scheme can be found in Rai (1987). 


5. BOUNDARY CONDITIONS 

The boundary conditions required when using multiple zones can be broadly classified 
into two types. The first type consists of the zonal conditions which are implemented at 
the interfaces of the computational meshes and the second type consists of the natural 
boundary conditions imposed on the surface and the outer boundaries of the 
computational mesh. The treatment of the zonal boundaries can be found in Rai 
(1986). The natural boundary conditions used in this study are discussed below. In 
particular radiating boundary conditions for the inlet and exit boundaries are 
presented. 


5.1. Airfoil Surface Boundary 

The boundary conditions on the airfoil surfaces are the “no-slip” condition and 
adiabatic wall conditions. It should be noted that in the case of the rotor airfoil, 
“no-slip” does not imply zero absolute velocity at the surface of the airfoil, but rather, 
zero relative velocity. In addition to the “no-slip” condition, the derivative of pressure 
in the direction normal to the wall surface is set to zero. 


5.2. Exit Boundary 

Two types of boundary conditions were used at this boundary. The first was a reflective 
boundary condition where the exit pressure was specified and three quantities are 
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extrapolated from the interior. The boundary conditions are 


P static constant. 

(13a) 

«i-o. ^=0 

(13b, c) 

ax ax 

II 

o 

(13d) 

ax 

Riemann variable, S 

- p/p y is the entropy and v is the 


and 


where R x - u + 2 cl(y — 1) is 
transverse velocity. This type of boundary condition reflects the pressure waves that 
reach the boundary back into the system. Two types of radiating boundary conditions 
were also implemented. The first was a one-dimensional boundary condition and its 
formulation (Bayliss & Turkei 1982) is shown below. It is assumed that, at the 
down-stream boundary, the flow is linear, that is, the unsteadiness can be considered a 
linear perturbation to a steady flow. This steady flow need not be axial. However, we 
can always rotate the coordinate system such that the x-axis is aligned with the 
underlying steady flow. In the rest of the analysis, it is assumed that the coordinates are 
so aligned. The linearized Euler equations far downstream are given by 
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= 0, (14) 


where u e , p e and c e are the underlying steady state velocity, density and speed of sound 
at the exit, whereas, p' , p' , u’ and v' are the perturbations in the pressure, the density, 
the velocity in the direction of the mean flow, and the velocity normal to the mean 
flow, respectively. Using the linearized x-momentum, y-momentum and energy 
equations, we can obtain 

p'"-2u e p' xt - {c 2 e - U 2 e )p' xx - c 2 e p' yy = 0. (15) 

Introducing the change in variables 


f ?■ v y 

and 

r = c e Vl - M\t + M£, 

equation (15) is transformed to 

P'rt-Pes-Pvi 7 = 0 - ( 16 ) 

Equation (16) admits solutions of the form 

p' —fix ~ £ cos 9 + Tj sin 9), (17) 

where 9 is the angle between the underlying mean flow and the axial direction 
( e ~ tan -1 (u«/i/4; u e = ^ul + vl). These solutions are planar waves propagating in the 
axial direction. Since in the present calculation the airfoils extend from p = -<» to + 
waves of this form are expected to exist. It should be noted that a boundary condition 
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based on equation (17) will be non-reflective for one-dimensional cases. Equation (1' 
suggests a differential operator of the form /£ — (B/3 t) + cos 0(3/9£) - sin 0(d/3rj 
which annihilates the functional form in equation (17). Hence the radiating boundai 
condition used is E£p — 0. which when transformed into the actual physical non-rotate 
coordinates give, for M^« 1, 


p't = Po=c co (l + M x ^){u' t + «*«(). (If 

Implementation of equation (18) is done by first setting u' x equal to zero. This 
consistent with the zeroth order extrapolation of the velocities at the boundary. Tb 
velocities u and v and the entropy p/p 7 are extrapolated from the interior. Th 
pressure is updated on the exit boundary by using equation (18). In this equation, th 
terms p^, and are obtained by circumferentially averaging p, u and v at the ex: 
at the previous time step. The exit sonic speed, c Xf is evaluated by using = VypZTpl 
where the value of p x is fixed and is equal to the exit pressure value used in th 
reflective boundary condition procedure. 

The boundary condition described above is essentially a one-dimensional boundar 
condition. Two-dimensional boundary conditions as developed by Giles (1988b) wer 
also implemented. As in the previous boundary condition, the flow at the exit : 
assumed to have small perturbations and hence be linearizable about an underlyin 
mean flow. The exit boundary conditions in terms of one-dimensional characteristi 
variables are 


3C 4 r r. , 3 , 

— - + {0 0 v x }~~ i 

dt By 
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,c 3 , 


(19a 


(19b 


where the transformation between the characteristic variables and the perturbatio 
variables is given by 
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Implementation of this boundary condition requires a knowledge of the underlying exit 
flow variables, p«, u v v , and p». The first three quantities are time lagged whereas 
the exit pressure, p*, is kept constant. 


5.3. Inlet Boundary 

Here again two types of boundary conditions were used. The first was the reflective 
boundary condition procedure wherein three quantities have to be specified. The three 
chosen are the Riemann invariant 


Ri = u + 


2c 

7-1 


U- = 


2 c_ K 

+ , 

7-1 


the total pressure 


P total 



1 + 




and the inlet flow angle, which in this case is equivalent to 

^iniet O' 


(21a) 


(21b) 


(21c) 


The fourth quantity needed to update the points on this boundary is also a Riemann 
invariant but is extrapolated from the, interior and is given by 



(21d) 


where 


#2 = w ~- 

7 - 1 

In the above equations the quantities u and v are the velocities in the x and y 
directions, p is the pressure and c is the local speed of sound. Specifying the total 
pressure at the inlet results in the boundary condition being reflective. 

A non-reflective or radiative one-dimensional boundary condition (Bavliss & Turkel 
1982) was also implemented. As in the case of the exit boundary, it is assumed that at 
the upstream boundary the flow is linear, that is. the unsteadiness can be considered a 
linear perturbation to a steady flow. Using an analysis that is very similar to that used 
in developing the radiating boundary condition for the exit boundary we obtain the 
following condition at the inlet: 


p/ = p-*(«-» - £-*)(«/ + u _*«;). 


( 22 ) 


At the upstream boundary, the Riemann invariant, R u and the flow angle, u in)ef = 0, are 
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still used. However, p total is replaced by the condition that at the inlet, the flow is 
isentropic, which gives 

p/p 7 = constant = p „Jp Z„, (23) 

and equation (21d) is replaced by equation (22). 

The radiating boundary condition described above is basically one-dimensional in 
nature. Two-dimensional inlet boundary conditions (Giles 1988b) were also imple- 
mented. The inlet boundary conditions in terms of the one-dimensional characteristic 
variables are 
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— < 
dt 


Cl 

c 2 > + 
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v 0 0 0 

0 V j(c + u) \(c-u) 
LO ?(c - u) v 0 
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dy 


Cl 

c 2 

c 3 

lc 4J 


= 0, 


dc 4 

dx 


= 0 . 


(24a) 


(24b) 


The transformation between the one-dimensional characteristic variables and the 
perturbation flow quantities are given by equation (20), with the quantities ( ) x 
replaced by the quantities ( )_*.. 


5.4. Upper and Lower Boundaries 

The computations reported in this study assume that the flow is spatially periodic in the 
y-direction. The spatial interval of .periodicity depends upon the airfoil count. (For 
example, in the 3-stator/4-rotor case, periodicity is imposed after every three stator 
airfoils and four rotor airfoils.) Further details regarding this boundary condition can 
be found in Rai (1987). 


6. RESULTS 

In this section, results obtained for the single-rotor/single-stator and four-rotor/three- 
stator configurations are presented. In particular, a comparison between the long-grid 
and short-grid computation with reflective boundary conditions and non-reflective 
boundary conditions will be made. In addition, the spectrum of the turbine tones and 
the variation of the amplitudes of the different modes in the far field will be presented 
for different cases. 

The dependent variables are non-dimensionalized with respect to the far upstream 
pressure (/?_,=) and density (p-*). The free-stream velocities are 

u -* - v_ x = 0, 

where = 0-07 is the inlet Mach number. The pressure ratio across the turbine 
(p static exil /p iotaii„ let ) 0'963. The rotor velocity was obtained so as to match the 
experimental flow coefficient (ratio of average inlet velocity to rotor speed) of 0*78 as 
given in Dring et al. (1982). The Reynolds number is 100,000/in. The kinematic 
viscosity was calculated using Sutherland’s law and the turbulent eddy viscosity was 
calculated using the Baldwin— Lomax model. The calculations were performed at a 
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constant time-step value of about 0-16 (this translates into 500 time steps for the rotor 
to move through a distance equal to the distance between two successive blades). 

6.1. Airfoil Surface Pressure Amplitudes 

The first comparison is made between the long and the short grid computations for the 
single-rotor/single-stator case. Reflective boundary conditions are used in both 
computations. Figure 2 shows the pressure amplitudes on the stator for the two cases. 
The symbols in this figure are the experimental data of Dring et al. (1982). The 
pressure amplitude (C p ) is defined as 

P max P min 

~ i T~ 

2PinIet^ 

where o) is the rotor velocity and p max and p min are the maximum and minimum 
pressures that occur over a cycle. A cycle corresponds to the rotor moving by a 
distance equal to the distance between adjacent rotor or stator airfoils. The pressure 
amplitudes obtained in the short-grid computation show most of the qualitative 
features that are found in the experimental results. However, the numerical data show 
a wider large-amplitude region than that found experimentally. In addition, the 
predicted peak is to the left of the experimental peak, and the pressure amplitude 
minimum on the suction side seen in the experimental results (x ~ -2-4) is absent in 
the computed results. The long grid computation yields an amplitude distribution closer 
to the experimental data. The position of the peak and the extent of the large 
amplitude region agree well with the data. 

The improvement obtained using the long grid is due to the large grid spacing in the 
far field region of the outer grid, which attenuates propagating modes. However, one 
penalty incurred in using this approach is the excessive computer time needed to obtain 
a periodic state. A typical short grid computation for a 3-stator/4-rotor case to 
converge to a time periodic state (including convergence in the tonal acoustics) is about 
20 cpu hours on a single processor of a CRAY YMP supercomputer. For a 
corresponding long grid computation, approximately five times as much computing 
time is required. Additionally, the time for convergence varies linearly with the number 
of airfoils, provided the extent of the upstream and downstream regions of the grid is 
the same. It should be mentioned that the time for convergence for flow quantities such 




Axial distance along stator surface 

Figure 2. Pressure amplitude on stator surface for (a) short grid, (b) long grid for a 1-stator/I-rotor case: 

O, suction surface; ♦, pressure surface. 
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Figure 3. Pressure amplitude on stator surface for a 1-stator/l-rotor case (using nonrefiective boundar 
conditions): O, suction surface; ♦, pressure surface. 


as pressure amplitudes on the airfoil surfaces and near field acoustics was considerabl 
less. 

Figure 3 shows the computed surface pressure amplitude distribution obtained usin; 
the short grid in conjunction with the one-dimensional non-reflective inlet and exi 
boundary conditions [see equations (18, 22)]. The agreement with the experimenta 
data is slightly better than that obtained on the long grid with reflective boundar 
conditions. It was found that the level of repeatability (solution periodicity in time 
with these boundary conditions was much better and the solution converged to a tim< 
periodic state faster. In addition, turbine operating conditions were maintained unlike 
in Rai (1990), where the use of non-reflective boundary conditions required an iterative 
process in which the Riemann invariant, R 2 , (specified at the exit) had to be variec 
until the proper average exit pressure was obtained. The flow coefficient and the 
turbine pressure ratio differed from that obtained using the reflective boundar; 
condition by less than 1%. (This will be shown later in Table 3.) Figure 4 shows the 



Axial distance along rotor surface 

Figure 4. Pressure amplitude on rotor surface with a (a) reflective boundary conditions and (b) 
non-reflective boundary conditions for a 1-stator/l-rotor case: O, suction surface; ♦, pressure surface. 
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Figure 5. Pressure amplitude on stator surface for 3-stator/4-rotor case with (a) a short grid and (b) a long 

grid: O, suction surface; ♦, pressure surface. 


pressure amplitudes on the rotor for a short grid with and without reflective boundary 
conditions. It is seen that the amplitudes obtained using the non-reflective boundary 
conditions are generally lower. 

In contrast to the single-stator/single-rotor case, the pressure amplitudes for the 
3-stator/ 4-rotor case did not differ much for the long or short grid or for the reflective 
or non-reflective boundary conditions. For the single-stator/single-rotor case, an 
acoustic analysis (Tyler & Sofrin 1970) shows that every harmonic could have 
propagating modes whereas the 3-stator/4-rotor case does not have any propagating 
modes for the first two blade passing harmonics. The first two harmonics do have 
decaying modes. In the original experimental configuration, there are 22 stator airfoils 
and 28 rotor airfoils. For this case also, the acoustic analysis does not predict any 
propagating modes for the first two blade passing harmonics. Since the higher 
harmonics are usually much smaller in magnitude, the unsteady pressures that reach 
the computational boundaries in the 3-stator/4-rotor case are much smaller than that 
for the single-stator/single-rotor case. Hence, it is expected that the reflective 
properties of the boundary conditions would play a smaller role in determining the 
unsteady pressures on the airfoils for the 3-stator/4-rotor case. Figure 5 shows the 
pressure amplitudes on the stator surface for the short and long grids with reflective 
boundary conditions. Figure 6 shows the pressure amplitude for the short grid with 



Figure 6. Pressure amplitude on stator surface for 3-stator/4-rotor case (using non-reflective boundary 
conditions): O, suction surface; ♦, pressure surface. 
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non-reflective boundary conditions. Clearly there is an improvement over the resul 
depicted in Figure 2(a) and a slight improvement over that depicted in Figures 2(b) ar 
3. The extent of the large amplitude region and the location of the pressure pes 
matches the experimental data. The slight improvement is due to a closer similarity 
the geometry with that of the experimental geometry. However, it should not l 
concluded that reflective properties of the computational boundaries are unimportai 
for the 3-stator/4-rotor case; they can still significantly alter the tonal acoustics in t! 
linear region of the flow. The pressure amplitudes on the rotor surface for th 
3-stator/4~rotor case also did not depend on the type of grid or the boundai 
conditions. The amplitudes on the rotor surface were very similar to that reporte 
previously by Rai & Madavan (1990). 

Besides the pressure amplitudes, phase information can also be obtained. The pha: 
of the low pressure peak on the stator suction surface [see Rai and Madavan (1990) f< 
details] for the 3-stator/4-rotor case did not depend on the type of grid or the boundai 
conditions. The numerical results compared well with experimental data of Dring er c 
(1982) and were similar to that reported by Rai & Madavan (1990). 


6.2. Far Field Linear Behavior 

The spectrum of turbine tones obtained from the computations is presented in th 
section. Recall that the Fourier modes predicted by the kinematical analysis ai 
denoted by p mn [equation (3)], where m and n are related as given by equation (5). Th 
values of a mn can be obtained by performing a Fourier decomposition of the pressui 
variation upstream and downstream of the turbine. The upstream resuits wei 
calculated at two chord-lengths upstream of the leading edge of the stator airfoils ar 
the downstream results were calculated at two chord-lengths downstream of the trailir 
edge of the rotor airfoils. Figure 7(a, b) shows the contribution of the axially aligne 
planar waves (m = 0) for the single-stator/single-rotor case. The ar-axis corresponds I 


xlCH xlCH 



Figure 7. Spectrum of the m - 0 mode (a) upstream of the stator and (b) downstream of the rot 

(1-stator/l-rotor). 
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XlO-5 X1Q-6 



Figure 8. Spectrum of the (a) m = 1 and (b) m = -1 modes downstream of the rotor (1-stator/l-rotor). 


harmonics of blade passing frequency and the y-axis corresponds to the computed 
coefficients, a mn . It is seen that, in general, the amplitudes of the higher harmonics are 
smaller than the amplitudes of the lower harmonics. Equation (6) predicts that all these 
harmonics propagate without decay. The numerical results conform with this prediction 
[see Rangwalla & Rai (1990)]. Note that the contribution due to the subharmonics of 
blade passing frequency is very small (by two orders of magnitude) compared to the 
harmonics of that frequency, thus leading to the conclusion that, for this mode, the 
kinematical interactions dominate. It was also found that the contribution of the 
non-planar modes ( m ¥= 0) upstream of the stator was less by at least an order of 
magnitude when compared to the planar modes. Figure 8 shows the m = 1 and m = -1 
modes downstream of the rotor. These modes are about an order of magnitude smaller 
than the pianar mode. The subharmonic content is very small as in the m — 0 case. 

The situation in the 3-stator/4-rotor case is different. For the planar case (m = 0), 
equation (5) predicts the existence of only the n — 3, 6, 9, . . . harmonics of blade 
passing frequency. Figure 9(a) shows the contribution of these planar waves upstream 
of the stator. Although the subharmonic content is more than in the single -stator/ 
single-rotor case, most of the energy is seen to lie in the n = 3 and n = 6 harmonics. 
Figure 9(b) shows the contribution of the planar waves downstream of the rotor. It 
should be noted that the pressure variations downstream of the rotor are measured in 
the rotor frame of reference. Hence for the m = 0 modes, the kinematical analysis 
predicts the existence of only the n = 4, 8, 12, . . . harmonics of blade passing frequency. 
Once again, we notice a low level of subharmonic noise, thus leading to the conclusion 
that for this case the kinematical interactions dominate. It should be noted that, for the 
m = 0 modes, the subharmonic noise of the single-stator/single-rotor case is of the 
same order of magnitude as the 3-stator/4- rotor case than that in the single- 
stator/single-rotor case. 

The m- 1 and m = — 1 modes upstream of the stator are shown in Figure 10(a, b). 
Substituting R — 4 and S = 3 in equation (5), the positive integer values that are 
possible for n when m = — 1 are n = 1, 4, 7, . . . and when m = 1, the positive integer 
values that n can take are n- 2, 5, 8, . . . We observe that the dominant frequencies 
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Figure 9. Spectrum of the (a) m = 0 mode upstream of the stator and (fa) the m — 0 mode downstream o 

the rotor for the 3-stator/4-rotor case. 


conform with the kinematical analysis. The same is true downstream of the rotor. Since 
the pressure variations downstream of the rotor are measured in the rotor frame o 
reference, R and S in equation (5) should be interchanged. Hence, for the m — 1 mode 
downstream of the rotor-stator pair, equation (5) predicts the existence o 
n - 3, 7, 11, . . . blade passing harmonics. Similarly, for the m - — 1 mode, equation (5' 
predicts that the blade passing , harmonics that can be present are given b) 

n = 1, 5, 9 The m — 1 and m ='— 1 modes downstream of the rotor are shown ir 

Figure 10(c, d). Once again we observe that the dominant frequencies conform with the 
kinematical analysis. 

The computations can also be used to study the propagation or decay of the various 
modes. To do this, either a long grid has to be used and the amplitudes of each mode 
calculated in the region of the grid where the solution has not suffered from numerica 
dissipation (due to grid coarseness) or a short grid with non-reflective boundary 
conditions should be used. A study of the numerical propagation or decay of the 
different modes has the advantage of determining the grid spacing required in the far 
field to maintain a propagating mode or to capture accurately the decay rate of £ 
decaying mode. Additionally, the effect of boundary conditions on the different modes 
can be assessed. 

Figure 11 (a, b) shows the effect of grid coarsening on propagating waves for the 
3-stator/4-rotor case. Figure 11(a) shows the amplitudes of three propagating modes 
[a Q3 , a 0 _ 6 and a Qt9 as given by equation (8)] upstream of the rotor-stator pair. As 
propagating modes, these amplitudes should remain constant. However, the amplitudes 
do decay as the grid spacing in the axial direction increases; (the symbols indicate the 
axial location of the grid points). As expected, the higher harmonics decay faster 
because they have smaller wavelengths. In the present case, the wavelength of the 
(0, 3) mode is approximately 21 T1 in. (536-2 mm), the wavelength of the (0, 6) mode is 
approximately 10-55 in. and the wavelength of the (0,9) mode is about 7-04 in. From 
Figure 11(a) we see that numerical dissipation sets in when there are five or fewer mesh 
points within a wavelength. Figure 11(b) shows the amplitudes of the propagating 
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Figure 10. Spectrum of the (a) m — 1 and (b) m — — 1 modes upstream of the stator and of the (c) m = 1 
and (d) m — — 1 modes downstream of the rotor (3-stator/4-rotor). 

modes (tf 0 . 4 and a 0 ,s) downstream of the rotor-stator pair. Unlike the upstream results, 
these amplitudes exhibit rapid variations near the rotor-stator pair. These oscillations 
eventually subside and the amplitudes remain constant till they monotonically decay 
because of grid coarsening. The rapid axial variation of the amplitudes immediately 
downstream of the rotor-stator pair is due to the nonuniformity of the mean flow. This 
nonuniformity is largely due to the velocity defects in the wakes of the rotor airfoils. 
The axial range over which the effect of this nonuniformity is felt depends upon the 
mode. It will be seen later that the rate of decay of the decaying modes and the axial 
wavelength of propagating modes can be significantly affected by the nonuniformity of 
the underlying mean flow. The effect of numerical dissipation due to grid coarseness 
downstream of the rotor-stator pair is similar to that observed in the upstream region. 
In the present case, the wavelengths of the (0, 4) and (0, 8) modes are approximately 
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Axial distance upstream of turbine 


Figure 11. Axial variation of the amplitudes of some propagating modes (a) upstream of the stator and (b) 
downstream of the rotor (3-stator/4-rotor); •. axial location of grid points. 


20-96 and 10-48 in., respectively. In Figure 11(b) it is seen that numerical dissipation 
affects the propagating modes when the number of mesh points within a wavelength 
are five or less. 

Figure 12(a, b) shows the instantaneous variation in the axial direction of the (0, 3) 
and (0, 4) modes. In the figure, the effect of grid coarsening on the waveform can be 
seen. Grid coarsening can affect the amplitude as well as the wavelength of the mode. 
The present numerical results seem to indicate that the wavelength variation due to 
grid coarseness is slight, as long as there are more than five grid points per wavelength. 


xlO-5 *10-5 




Figure 12. Instantaneous axial variation of the (a) (0,3) mode upstream of the stator and (b) (0,4) mode 
downstream of the rotor (3-stator/4-rotor); •, axial location of grid points. 
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Table 1 


Axial wavelengths of some propagating modes (upstream); 3-stator/4-rotor case. 

(Theory in parentheses) 



-2 

-1 

0 

1 

2 

n = 1 
n—2 
n~ 3 

decaying 

decaying 

20-54(21-11) 

decaying 

decaying 

n =4 


20*81(21*63) 



decaying 

n = 5 

decaying 



14*96(15*09) 


n = 6 
n= 1 


9*46(9*81) 

10*49(10*55) 


14-15(14*48) 

n = 8 
n= 9 
n = 10 

10*57(10-82) 

6*35(6*58) 

7*03(7*04) 

8*11(8*41) 

7*16(7*55) 


However the amplitude can be significantly affected by grid coarseness and decays to 
about half its value when there are six grid points per wavelength. 

Comparisons between the theoretical axial wavelengths and that obtained numeri- 
cally are shown in the Tables 1 and 2 for the 3-stator/4-rotor case. 

Table 1 shows comparisons upstream of the rotor-stator pair. In general, the results 
are good. The differences are less than the maximum grid spacing in the region where 
the wavelengths were measured. Table 2 shows a similar comparison downstream of 
the rotor-stator pair. The numerical results are generally in fair agreement with the 
theoretical predictions for all (m, n) modes where m < 0. The differences are of the 
same order as the maximum grid spacing in the region where the wavelengths were 
measured. However, the calculated wavelengths of the im , n ) modes where m >0 do 
not agree well with theoretical predictions. It should be recalled that the theoretical 
prediction of axial wavelengths [equation (8)] was obtained under the assumption of a 
uniform mean flow. However, downstream of the rotor-stator pair, the underlying 
mean flow is not uniform due to the wakes of the airfoils. It is believed that this 
nonuniformity in the mean flow is the main reason for the discrepancy between the 
numerical results and theoretical predictions. 

One objective of the present study is to see if non-reflective boundary conditions can 
be used along with a short grid to predict the tonal acoustics present in rotor-stator 

Table 2. 


Axial wavelengths of some propagating modes (downstream); 3-stator/4-rotor case. 

(Theory in parentheses) 



-2 

-1 

0 

1 

2 

n = 1 
n— 2 
n = 3 
n =4 

decaying 

decaying 

23*51(20*96) 

decaying 

decaying 

n = 5 


23*78(26*91) 




n = 6 

decaying 




decaying 

n=l 
n= 8 
n= 9 
n = 10 

11*89(13*46) 

10*81(10*47) 

11*65(10*48) 

18*38(13*39) 

20*81(11*12) 
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Figure 13. Comparison between a long grid solution and a short grid solution with reflective boundai 
conditions for (a) the (0, 3) mode upstream and (0, 4) mode downstream of the rotor-stator pair, and (b) th 

(1,2) mode upstream and (1,3) mode downstream of the rotor-stator pair. , . Long grid, 

, short grid. 


interactions. Three different boundary conditions were compared. The first boundar 
condition was a reflective boundary condition as given in Rai & Madavan (1990). Thes 
boundary conditions have been widely used in the numerical simulations of rotor 
stator interactions. However, they are not adequate in the study of tonal acoustic 
because of their reflective properties. Figure 13(a, b) shows the comparison of shoi 
and long grid calculations with reflective boundary conditions for the 3-stator/4-rotc 
case. In Figure 13(a) the axial variation of the amplitudes of two propagating mode 
are shown. The effect of using reflective boundary conditions can be seen in this figuri 
At the exit boundary of the short grid, the amplitudes become zero whereas at the ink 
boundary, the amplitudes are very small. 

The short and long grid solutions also do not compare well within the flow domaii 
The short grid solutions do not show an axial region where the amplitudes remai 
constant. This is because a reflective boundary condition reflects any propagating moc 
back into the flow domain. As expected the long grid solution does exhibit an axi; 
region over which the amplitudes remain constant. The reason for this is that tl 
coarseness in the grids near the inlet and exit boundaries essentially dissipates tl 
propagating waves, thus minimizing the effects of reflection. In contrast to tl 
propagating modes, the effect of the reflective boundary conditions on the decayir 
modes is only significant near the boundaries. Figure 13(b) shows the axial variation < 
the (1,2) mode upstream of the stage and the (1,3) mode downstream of the stag 
The results show that the amplitudes remain unaffected in the near field region of tl 
airfoils. It should be recalled that these results are for the 3-stator/4-rotor case. For th 
case the boundary conditions did not significantly affect the pressure amplitudes on tl 
airfoil surfaces. The amplitudes on the airfoil surfaces are mainly composed of tl 
lower harmonics. For the 3-stator/4-rotor case, the lower harmonics decay wii 
increasing axial distance from the airfoils. Reflective boundary conditions do refle 
these modes but the effect of the reflection is confined to the region near tl 
boundaries. 
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Table 3. 

Converged operating conditions 



Experimental 

Reflective B.C. 

1-D B.C. 

2-D B.C. 

Pstatic a „ii 
P total in | c , 

V elocity i n i e t 

P totali nle , 

0-963 

0- 083 

1- 0034 

0-963 

0- 084 

1- 0034 

0-9636 

0- 083 

1- 0031 

0-963 

0- 084 

1- 0034 


modes for the short grid case are stmilar to that for the 

levels are a bit different. This difference short and long 

operating conditions as shown in Table . rnnd ; tions a re larger than that obtained 
grid solutions with one-dimensiona oun ary , a i s0 show some reflectivity 

with the two-dimensional ^*“"" 0 near the upstream 
at the boundaries, as evidenced by variation of the (1 2) mode and 

boundaries Figure 15(a,b) » the propagating 

the downstream axial variation of the (1, 3) m nea r-field region of the 

modes, the amplitudes of these modes remain unacted. he near orid 8 there are 

airfoils. At the inlet and exit computawna boundaries of ^shon ^ ^ 

£h*s~i sjzss z 

theTame modes remain nearly unaffected 

in the near-field regions of the airfoils. „vnnnentiallv in the upstream and 

The amplitudes Of the defying . modes r ^ o “^ dep i ds upon the underlying 

downstream directions [equation (11}]- f , mode as ojven in equation 

mean flow and the temporal and spat, a ^enc^of ttemode asjven^ ^ ^ 

perturbations ^o^a s^dy uniform mean flow], 

VOX ssssss <s 2: - 

upstream and downstream of the rotor-stator pain It shou d be note q 

(ID and (12) yield only the rate of decay and * s ^ U of the e mSdes 


Table 4, 

Mean flow quantities 



Upstream 

Downstream 

Axial velocity, U x 
Transverse velocity, U, 
Sonic velocity, c 

0- 0842535 
0-00 

1- 1829303 

0-0864004 

0- 0701666 

1- 1774635 
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Figure 16. Comparison between numerical and theoretical decay rates upstream for (a) mode (1,2) and 
(b) mode (-1,1) upstream of the stator and for (c) mode (—1,1) and (d) mode (1,3) downstream of the 

rotor: , numerical; , theoretical. 

solution deviates from the theoretical exponential decay because of the very coarse grid 
in this region. 

Downstream of the rotor-stator pair, a difference between the numerical results and 
the linear theory is seen [Figure 16(c, d)]. The axial variation of the natural logarithm 
of amplitude of the ( — 1, 1) mode (log a_ ia ) is shown in Figure 16(c). Even though the 
overall decay matches that of the theoretical exponential decay, the numerical variation 
of the amplitude is not a “pure” exponential. Rather, the axial variation of the 
computed amplitude is exponential as well as oscillatory in nature. This behavior was 
observed for all (/, k) modes, where l < 0. Figure 16(d) shows the axial variation of the 
natural logarithm of the amplitude of the (1,3) mode. The decay of this mode is 
exponential, however, the decay rate does not match the theoretical prediction. This 
difference between the numerical and theoretical decay rates was observed for all (l, k) 
modes where / > 0. It should be recalled that one of the assumptions underlying the 
theoretical results is a uniform mean flow. However, downstream of the rotor-stator 
pair, the underlying mean flow does deviate from a uniform flow because of the wakes 
of the rotor airfoils. This deviation is much more than the deviation upstream of the 
rotor-stator pair due to potential effects. It is believed that the nonuniformity in the 
mean flow is the main reason for the discrepancy between the numerical results and 
theory. 


7. SUMMARY 

This study focuses on the numerical computation of tones in rotor-stator interactions, 
the numerical predictions are obtained by solving the thin-layer Navier-Stokes 
equations on a system of patched grids. The mode-cotent of interaction tone noise is 
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obtained for two different airfoil counts and is shown to conform with a kinematica. 
analysis of the flow. In addition, the propagation characteristics of different modes are 
compared with the predictions of linear theory. Numerically computed pressure 
amplitudes on the surface of the airfoils are compared with experimental data. 

The effects of both reflective and non-reflective boundary conditions on the 
calculated flow field were assessed. It was found that for the short-grid calculations 
non-reflective boundary conditions had to be used in order to predict the tona 
acoustics in the flow field. Use of non-reflective boundary conditions becomes more 
important for those cases where there is a high energy content in the propagating 
modes. (The single-stator/singie-rotor case for example, had propagating lowei 
harmonics. These harmonics had high energy content. Reflection of these harmonic: 
from the computational boundaries resulted in a degradation of the computed pressure 
amplitudes on the airfoil surfaces.) It was also shown that reflective boundary 
conditions can be made to behave essentially as non-reflective boundary condition: 
with the use of appropriately long grids in the upstream and downstream regions. The 
long grid approach in conjunction with an increasing grid cell size in the far upstrean 
and downstream regions results in numerical dissipation of propagating modes, thu: 
avoiding the problems due to reflections at the computational boundaries. 
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Abstract 

This paper deals with the application of a three- 
dimensional, time-accurate Navier-Stokes code for pre- 
dicting the unsteady flow in an advanced transonic 
turbine. For such advanced designs, prior work in 
two dimensions has indicated that unsteady interac- 
tions can play a significant role in turbine performance. 
These interactions affect not only the stage efficiency 
but can substantially alter the time-averaged features 
of the flow. This work is a natural extension of the 
work done in two dimensions and addresses some of 
the issues raised therein. These computations are be- 
ing performed as an integral part of an actual de- 
sign process and demonstrate the value of unsteady 
rotor-stator interaction calculations in the design of 
turbomachines. Results in the form of time-averaged 
pressures and pressure amplitudes on the airfoil sur- 
faces are shown. In addition, instantaneous contours 
of pressure and Mach number are presented in order to 
provide a greater understanding of the inviscid as well 
as the viscous aspects of the flowfield. Relevant sec- 
ondary flow features such as cross-plane contours of to- 
tal pressure and span-wise variation of mass-averaged 
quantities are also shown. 

Introduction 

The traditional design of new turbines has re- 
lied upon empirical correlations, extensive experimen- 
tal data, and a technology data base comprising pre- 
vious designs. This design process has proven very 
reliable for new designs that do not deviate very much 
from those in the existing data base. However, a more 
general predictive capability is needed when the oper- 
ating conditions of a new design demand radical devi- 
ations from the data base. 

Considerable progress has been made in using 
computational fluid dynamics (CFD) to predict flows 
within turbomachines. Much of the early work has fo- 
cused on predicting the flow in airfoil cascades. An 
extensive body of experimental and numerical results 
in the literature deals with a wide variety of two-and 
three-dimensional cascade geometries. While such meth- 
ods of analysis of flows in isolated airfoil rows have 
helped improve our understanding of flow phenomena 
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in turbomachinery and have gained widespread accep- 
tance in the industrial community as a design tool, 
they do not yield any information regarding the un- 
steady effects arising out of rotor-stator aerodynamic 
interaction. However, it is becoming increasingiy im- 
portant to consider interaction effects in the design of 
new generation turbines. This has come about due to 
the constraints of low weight, small size, high specific 
work per stage, high efficiency, and durability, which 
results in very high turning angles and unconventional 
airfoil shapes, potentially giving rise to nonlinear un- 
steady interactions. 

In the past few years, advanced transonic tur- 
bines have been designed by Pratt and Whitney in 
support of the Consortium for CFD Application in 
Propulsion Technology sponsored by NASA Marshall 
Research Center. These turbines are characterized by 
very high flow turning angles (160° per stage) and rel- 
atively high loading coefficients. The current status 
of turbine design technology is shown in Fig. 1 (pri- 
vate communication, L. Griffin, NASA Marshall Space 
Flight Center). The figure shows the turning angles 
and turbine loading coefficients of some existing tur- 
bine designs. The figure shows two designs (the G 3 T 
and the G 2 OT ) that have turning angles of 160° which 
is 20° higher than the traditional design limit of 140°. 

Numerical methods that simulate the unsteady 
flow associated with rotor-stator configurations have 
been developed in recent years. References 1—3 present 
a zonal approach for solving the unsteady, thin-layer, 
Navier-Stokes equations for rotor-stator configurations 
in a time-accurate manner, both in two and three di- 
mensions. The present work is an application of the 
three-dimensional rotor-stator code described in Ref. 
3, to evaluate the design of the advanced Gas Genera- 
tor Oxidizer turbine designed by Pratt and Whitney, 
and will henceforth be referred to as the G 2 OT. It 
should be mentioned that the first application of an un- 
steady two-dimensional Navier-Stokes solver for design 
purposes was carried out for the G 3 T (Ref. 4). The nu- 
merical predictions were obtained at a constant radius 
corresponding to the midspan of the rotor airfoils. The 
primary issue was the effect of unsteady interactions 
on boundary layer separation. The results from the 
unsteady two-dimensional analysis led to design mod- 
ifications (Ref. 4) and provided the designers with a 
better understanding of the physics of the flow. The 
results also validated the concept of using high turn- 
ing angles and high specific work per stage. One out- 


come was the possibility of using single-stage turbines 
for certain applications, thus considerably simplifying 
the design process. This directly influenced the de- 
sign of the G 2 OT. The unsteady two-dimensional code 
was again used to aid in the design of the G 2 OT and 
was reported in Ref. 5. The flow was predicted at 
a constant radius which was equal to the midspan of 
the rotor airfoil. Results were obtained for two power 
settings (100% and 70%) and it was found that the 
turbine loads were within the tolerances specified by 
design requirements and were acceptable. The two- 
dimensional analysis used in Ref. 4-5 contained quasi- 
three-dimensional source terms to account for stream 
tube contraction effects (Ref. 6). The numerical algo- 
rithm was an extension of that previously reported in 
Ref. 7. 

One drawback of the two-dimensional analysis 
is that it is not complete. Since the flow is three- 
dimensional, the issue of secondary flow influencing 
flow features such as strength and position of the shocks 
has to be addressed. Other issues such as the effect of 
unsteady interactions on the end-wall boundary layers 
have to be assessed. Hence a three-dimensional inter- 
action study was initiated. 

In this paper, the three-dimensional as well as 
some two-dimensional results for the G 2 OT will be 
presented. The two-dimensional results were obtained 
for two power settings (100% and 70%) whereas the 
three-dimensional results were obtained for only the 
100% power setting. Comparisions will be made wher- 
ever possible. In particular, it was found that there 
were similarities as well as differences between the two- 
dimensional and three-dimensional results. The over- 
all loading on the airfoils obtained from the three- 
dimensional analysis at midspan compared fairly with 
the two-dimensional predictions. However, some of the 
details such as strength and positions of the shocks dif- 
fered. This also resulted in weaker unsteady interac- 
tion predictions by the three-dimensional calculations. 

Two grid systems (one with twice as many points 
in the radial direction than the other) were used for a 
limited grid independance study. Each grid system 
contains multiple patched and overlaid grids as de- 
scribed in Ref. 3. These grids can move relative to 
one another to allow for the relative motion of the ro- 
tor airfoils with respect to the stator airfoils. 

Numerical Method 

The numerical method solves the unsteady, three- 
dimensional, thin-layer Navier-Stokes equations. The 
Navier-Stokes equations in three dimensions are nondi- 
mensionalized and transformed to a curvilinear time- 
dependent coordinate system, and a thin-layer approx- 


imation is then made. The unsteady, thin-layer, Navier- 
Stokes equations are solved using an upwind-biased 
finite-difference algorithm. The method is third-order- 
accurate in space and second-order-accurate in time. 
Several iterations are performed at each time step, so 
that the fully implicit finite-difference equations are 
solved to ensure a time-accurate solution. Further de- 
tails of the method can be found in Ref. 3. 

Boundary Conditions 

The boundary conditions required when using 
multiple zones can be broadly classified into two types. 
The first are the zonal conditions which are imple- 
mented at the interfaces of the computational meshes, 
and the second are the natural boundary conditions 
imposed on the surface and the outer boundaries of 
the computational mesh. The treatment of the zonal 
boundary conditions can be found in Ref. 2. The nat- 
ural boundary conditions used in this study are dis- 
cussed below. 

Airfoil Surface Boundary 

The boundary conditions on the airfoil surfaces 
are the “no-slip” condition and adiabatic wall condi- 
tions. It should be noted that in the case of the rotor 
airfoil,' “no-slip” does not imply zero absolute velocity 
at the surface of the airfoil, but rather, zero relative 
velocity. In addition, the derivative of pressure in the 
direction normal to the wall surface is set to zero. 


Exit Boundary 


The flow in the axial direction is subsonic at the 
exit boundary and hence only one flow quantity has to 
be specified. The flow quantity chosen in this study 
is the exit static pressure as a function of radius. To 
completely specify the flow variables at the boundary, 
four other flow quantities are extrapolated from the 
interior. The four chosen are the Reimann invariant, 

2c 

Ri = u 4 - 

7 - 1 


the entropy, 



and the velocities in the transverse directions. One dis- 
advantage of this type of boundary condition is that 
the pressure waves that reach the boundary are re- 
flected back into the flow domain. However, this bound- 
ary condition was chosen since, in general, it provides 
greater control on the turbine operating conditions 
and results in the correct pressure drop and mass flow 



through the turbine. 


Geometry and Grid System 


Inlet Boundary 


The flow at the inlet boundary is subsonic. Four 
quantities need to be specified at this boundary. The 
four chosen were the Reimann invariant, 

2 c 

Ri ~ u H 7 

7-1 

the total pressure as a function of radius, 

7 — 1 y 

Ptotal — Pinlet (l 4* ^ ^ inlet ) 

and the inlet flow angles, 

Vinlet 



and 


Winlet 
n inlet 


— tan(<j}) 


The fifth quantity needed to update the points on this 
boundary is also a Reimann invariant that is extrapo- 
lated from the interior and is given by 


In the above equations, the quantities u and v and w 
are the velocities in the axial (x) tangential (9) and 
the radial (r) directions, p is the pressure and c is 
the local speed of sound. Specifying the total pressure 
at the inlet results in a reflective boundary condition, 
but together with the specification of the exit static 
pressure, has the advantage of determining uniquely 
the turbine operating conditions. 


Periodic Boundaries 


Turbomachines are designed with unequal airfoil 
counts in the stator and rotor rows in order to mini- 
mize vibration and noise. A complete viscous simula- 
tion including all of the airfoils in the stator and rotor 
rows is yet impractical in a design environment. The 
approach used here is to assume that the ratio of the 
number of stator to rotor airfoils is a ratio of two small 
integers. This is achieved by scaling the stator or the 
rotor geometries such that the blockage remains the 
same. Periodicity conditions are then imposed over 
the composite pitch. For the case of the G 2 OT tur- 
bine, the number of stator airfoils is 20 and the number 
of rotor airfoils is 42. By changing the number of sta- 
tor airfoils to 21 and rescaling the stator airfoils by a 
factor of 20/21, a stator to rotor airfoil count of 1 to 
2 is achieved. The calculation assumes that the flow 
exhibits spatial periodicity over one stator airfoil and 
two rotor airfoils. Note that the pitch of one rescaled 
stator airfoil is equal to the composite pitch of two 
rotor airfoils. 


A schematic diagram of the G 2 OT is shown in 
Fig. 2. This is a single stage turbine that is designed 
to operate in the transonic regime. It is characterized 
by very high turning angles and high specific work. 

Figures 3a-b show the system of overlaid grids 
used to discretize the flow domain. The figure shows 
the fine grid with 51 grid points in the span wise direc- 
tion. Figure 3a shows the grid at the midspan. Each 
airfoil has two zones associated with it; an inner zone 
and an outer zone. The inner zone contains an O-grid 
that is generated using an elliptic grid generator. This 
grid is clustered near the airfoil surface in order to re- 
solve the viscous effects. The outer zone is discretized 
with an H-grid and is generated algebraically. The in- 
ner and outer grids overlap one another. This position- 
ing of the inner and outer grids facilitates information 
transfer between the two zones. The outer H-grids of 
the stator airfoils and rotor airfoils overlap and slip 
past each other as the rotor airfoils move relative to 
the stator airfoils. Figure 3b shows the surface grid 
(minus the casing). Here, the grid in the tip clearance 
region can also be seen. This grid was also generated 
by means of an elliptic grid generator and maintains 
metric continuity with the inner O-grid. The fine grid 
contains approximately 940000 grid points whereas the 
coarse grid has half as many. 

Results 

It should be mentioned that the numerical method 
has been validated both in two and three-dimensional 
applications. In particular, the ability to predict the 
time-averaged pressures and pressure amplitudes on 
airfoil surfaces and total pressure losses in airfoil wakes 
have already been demonstrated for turbines as well as 
for compressors (see Refs. 1-4, 7, 8). 

G 2 OT Two-Dimensional Computations 

A brief description of the two-dimensional results 
(Ref. 5) will first be presented for the purpose of 
comparison with the three-dimensional results. Two- 
dimensional predictions were obtained for two power 
settings. The first setting is at 100% power and the 
second is at 70% power. The operating conditions for 


the two power settings are shown in Table I. 


G 2 OT Three-dimensional results 


Inlet Mach No. 

100% Power 
0.46 

70% Power 
0.54 

Inlet Reynolds No. 

2.6 x 10 6 /inch 

1.1 x 10 & /inch 

RPM 

7880 

6232 

Inlet Ptotai 

542.77psia 

313.82psza 

Exit Pstatic 

200.00psia 

144.5psia 

Inlet T to t a i 

1307. 02°R 

1080. 88°R 


Table I. Turbine operating conditions 


The results for the three-dimensional computa- 
tions of the G 2 OT for the 100% power setting are pre- 
sented in this section. These results were obtained 
by integrating the governing equations and boundary 
conditions described earlier. A modified version of the 
Baldwin-Lomax turbulence model (Refs. 9-11) was 
used to determine the eddy viscosity. The modifica- 
tion involves the use of a blending function that varies 
the eddy viscosity distribution smoothly between the 
blade and endwall surfaces. Further details can be 
found in Refs. 10-11. The kinematic viscosity was 
calculated using Sutherland’s law. 


Static Pressure Variation on Airfoils Fig- 
ures 4 and 5 show the time-averaged and unsteady en- 
velope of static pressure on the airfoil surfaces for the 
two different power settings, respectively. The pres- 
sure coefficient on the surface of the stator airfoils in 
this case is defined as 

C i _ Pstatic 

p ~ ~ 

Ptotali n ut 

where p sta tic is either the time-averaged static pres- 
sure on the surface of the airfoil (to obtain the time- 
averaged pressure distribution) or the maximum or 
minimum pressure over a cycle (which results in the 
pressure envelope). The time averaging is performed 
over a cycle which corresponds to the rotor airfoils 
moving through two airfoil pitches. The pressure co- 
efficient on the surface of the rotor airfoils is defined 
as 

Pstatic 

p — ~ 

Ptotal(relative) roi0 riniet 

Here, the pressure is normalized with respect to the 
time-averaged relative inlet total pressure to the ro- 
tor rows. The figures show that the results of the 
two power settings are qualitatively similar. The pre- 
dicted pressure amplitudes are slightly smaller for the 
70% power setting than for the design setting (100% 
power). The pressure distribution indicates a weak 
(nearly stationary) shock on the suction surface of the 
stator airfoil that impinges on the rotor suction surface 
near the leading edge (Ref. 5). It is this shock that 
accounts for the moderately high pressure amplitudes 
near the leading edge of the rotor airfoils. The pressure 
distributions also indicate a shock near the trailing 
edge of the rotor airfoils. This second shock is nearly 
stationary with respect to the moving rotor airfoils. It 
should be noted that these are two-dimensional results 
at constant radius. In the three-dimensional case, the 
interaction effects are found to be less severe due to 
the relaxation effects of the spanwise direction. 


Static Pressure Variation on Airfoils Fig- 
ures 6-8 show the time-averaged and unsteady enve- 
lope of static pressure on the stator and rotor airfoils 
at three spanwise locations. Figures 6a-8a show the 
pressure variations on the stator airfoil at the hub, the 
midspan and at the casing, whereas Figs. 6b-8b show 
the variations at the hub, the midspan and at the tip 
of the rotor airfoils. These results were obtained by 
the fine grid calculations and do not show any signif- 
icant differences when compared with those obtained 
from the coarser grid. The level of unsteadiness on 
the stator airfoils is small compared to that on the 
rotor airfoils. The amplitudes also are smaller at the 
casing than at the hub. The figures seem to indicate 
the existence of a weak shock (made clearer by contour 
plots) on the suction surface of the stator near the hub. 
The pressure amplitudes on the rotor airfoil are larger. 
The rotor airfoils are unloaded considerably at the tip. 
However, it was found that this is very localized near 
the tip region and is not very critical. The predicted 
pressure amplitudes of the three-dimensional results at 
midspan, are smaller than the two-dimensional results. 
This is mainly due to the difference in the strength of 
the predicted axial gap shock. 

Figures 9a-b show the comparisions of the time- 
averaged pressures between the three-dimensional and 
the two-dimensional results. On the stator airfoil, the 
two-dimensional calculations predict a lower unload- 
ing at the airfoil nose than that shown by the three- 
dimensional calculations. It should be noted that the 
two-dimensional calculations were performed on a sur- 
face of constant radius. Quasi-three-dimensional source 
terms associated with stream-tube contraction were 
included in the calculation, but the terms associated 
with radius variation were not. To properly account 
for these terms, the two-dimensional calculations would 
have to be performed on a cylindrical surface with an 
axially varying radius. The overall loading on the ro- 
tor airfoils compares better. However, the details are 
different. In particular, the position and strength of 



the trailing edge shock on the suction surface is differ- 
ent. Also, other details that were present in the two- 
dimensional calculations, such as a large local vari- 
ation on the suction surface, is absent in the three- 
dimensional results. 

Instantaneous Mach. Number Contours Fig- 
ures lOa-c show instantaneous Mach number contours 
at 20%, 50% (midspan), and 80% of span respectively. 
These results were obtained for the fine grid system. 
At 20% of span, the shock near the trailing edge of 
the rotor airfoil can be seen. At this location, an axial 
gap shock was also seen. However, unlike the two- 
dimensional calculations the three-dimensional calcu- 
lations predict an intermittent axial gap shock. At 
midspan (Fig. 10b), the shock in the axial gap region 
is much weaker, whereas there is no shock at the down- 
stream location. In fact, it was found that the radial 
extent of the axial gap shock varied with time with a 
maximum extent of about 50%. Figure 10c shows the 
Mach contours at 80%. The contours seem to indicate 
that there might be unsteady separation on the suction 
surface of the rotor airfoils. Recall that the turning an- 
gles in this turbine are very high, and there is a con- 
cern about massive boundary layer separation under 
the influence of unsteady interactions. The numeri- 
cal calculations do not predict massive boundary layer 
separation, as indicated by the Mach number contours, 
thus increasing the confidence in the design. 

Instantaneous Static Pressure Contours Fig- 
ures lla-c show the instantaneous pressure contours 
at 20% 50% (midspan) and 80% of span respectively. 
These contours basically highlight the inviscid features 
of the flow. As expected, the rotor shock near the hub 
can be seen. 

Mass- Averaged Quantities versus Span Fig- 
ure 12 shows the mass-averaged meridonial angle ver- 
sus normalized span at four different axial stations. 
The axial stations correspond to the inlet of the tur- 
bine, the midgap, half a chord downstream of the rotor 
airfoil, and about one and a half chord lengths down- 
stream of the rotor airfoil. The figure does show that, 
to a large extent, the flow turns about 160° through 
the stage, however, it also shows a region of under- 
turning at the midspan. Figure 13 shows the mass- 
averaged radial pitch angle. Recall from the schematic 
of the G 2 OT (Fig. 2) that the casing angle is -30° at 
the inlet, and is positive aft of the rotor (approximately 
11.75°). This is reflected in the mass averaged pitch. 
Figures 14-17 show the variation of the mass-averaged 
Mach number, the relative Mach number (relative with 
respect to the rotating rotor airfoils), the absolute to- 
tal pressure and the relative rotational total pressure. 


One surprising aspect of the results is the local increase 
in total pressure losses at the midspan. 

Time- Averaged Contours Figures 18a-b show 
the time-averaged contours of the relative rotational 
total pressure at midgap and half a chord length down- 
stream of the rotor airfoils. The circumferential extent 
of Fig. 18a equals the circumferential pitch between 
two successive stator airfoils whereas that of Fig. 18b 
equals the pitch between two rotor airfoils. Also, it 
should be noted that the time averaging is done in two 
different frames of reference. At the midgap (Fig. 18a) 
the frame of reference is stationary, whereas, down- 
stream of the rotor airfoils (Fig. 18b) it is rotating. 
The contours at the midgap do show the expected 
(nearly uniform in span) stator wake along with the 
hub and casing secondary flows. However, aft of the 
rotor blades, at the midspan, a region of slightly higher 
losses exists. This was also observed in the mass- 
averaged numerical data (Figs. 16-17). Figures 19a-b 
show the time- averaged contours of Mach number rel- 
ative to the rotor airfoils at the same axial location. 
The relative Mach number of the flow is subsonic at 
the midgap, but downstream of the rotor it becomes 
supersonic and eventually shocks. 

Summary 

A detailed numerical calculation of the three- 
dimensional unsteady flow in an advanced gas gener- 
ator turbine is presented. The computational results 
are obtained by solving the three-dimensional, thin- 
layer, Navier-Stokes equations on a system of overlaid 
grids. The numerical results do capture many aspects 
of the flow that could aid in the understanding of the 
flow. In addition, the results do not indicate any sig- 
nificant boundary layer separation, (an object of con- 
cern). The unsteady loadings were found to be within 
acceptable limits. 

The present results indicate that a proper un- 
derstanding of the unsteady interaction effects could 
play an important role in the design of advanced gas 
generator turbines. 
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Fig. 1 Current status of turbine design technology. 
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Fig. 2 Schematic of G 2 OT flow path. 



Fig. 3a Rotor-stator airfoil grid systems for the G 2 OT, 
at mid-span 



Fig. 3b Rotor-stator airfoil grid systems for the G 2 OT , 
surface grid. 
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'ig. 6 Static pressure variation on (a) stator and (b) rotor airfoil surfaces for the G 2 OT (Hub). 
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Fig. 7 Static pressure variation on (a) stator and (b) rotor airfoil surfaces for the G 2 OT (Mid-span). 
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Fig. 8 Static pressure variation on (a) stator (casing) and (b) rotor airfoil (tip) for the G 2 OT. 
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Fig. 9 Comparision between two- and three-dimensional predictions of time-averaged static pressures (a) 
stator and (b) rotor airfoil surfaces (Mid-span). 







-ireds %08 (o) ptre ireds %og (q) ‘reds ‘reds %08 (o) pre reds %og (q) ‘reds %qs 

%0Z O 2 ) W wnsswd jo sino^uoo snoauB^u-B^sui xi (*) P 2 Jsqumu qrej^ jo sinoijuoo enosirextrexstii QI '2ij 











1 ,0 _ 


0.9'IL* 

ft 

0 . 8- r 

| 0.7- f- 


— f. 


0.31- 


Stadoos 

__ Met to stator 

irfew 

— *- Further art 


Radial span 


Radial span 


Fig. 12 Mass-averaged meridonial angle as a function Fig. * 14 Mass-averaged absolute mach number 
of sDan. function of span. 
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Fig. 13 Mass-averaged pitching angle as a function of Fig. 15 Mass-averaged relative mach number as a func 
sDan tion of span. 
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Fig. 19a Time- averaged relative madb. number at mid- 
gap. 



Fig. 19b Time-averaged relative mach number aft of 
the rotor airfoil. 


